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Abstract 

Herein a genetic algorithm for optimising the design of layered 2 D het- 
erostructnre is proposed. As a proof-of-concept it is applied to Sb2Te3-GeTe 
phase-change material superlattices, and the resulting lowest energy struc¬ 
ture is grown experimentally. The similarity of the computational and ex¬ 
perimental structures is verified with the comparison of XRD spectra. The 
structure is found to be within 0.92 meV/at. from the energetically most 
favorable known structure for Ge2Sb2Te5. 

1 Introduction 

The chalcogenide phase-change materials (PCM) are widely used in data storage due to 
their unique nonlinear crystallization rate as a function of temperature, and the high 
optical and electrical contrast between the structural states[I]. These materials are 
stable in either phase for years at room temperature, but switch on the nanosecond 
timescale at elevated temperatures (typically 150-300 °C). The material properties, such 
as crystallization temperature, crystallization rate at elevated temperature, and data 
retention at room temperature are generally tuned by changing the composition of the 
material. However, changing the PCM composition changes all of the properties so that 
fine-tuning one property is likely to affect the others. Recently, it has been reported that 
it is possible to build layered superlattice structures [l2] of phase-change materials instead 
of simple alloys. This opens new degrees of freedom for optimization as the properties 
of these structures can be tuned by changing the layer thicknesses or sequence while 
keeping the composition constant [8]. 

Another class of layered materials is the van der Waals heterostructuresjlj, which 
comprise multiple layers of 2D crystals bound together by relatively weak van der Waals 
forces. Such materials include graphene, hexagonal boron nitride, and molybdenum 
disulfide. Typically these structures are made by exfoliating the 2D crystals layer by 
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layer, and stacking the layers into the heterostructure one by one. This is currently an 
active field of research since many material properties are expected to differ in atomically 
thin materials when compared to bulk. Combining 2D crystals with novel properties to 
superlattices can lead to vdW materials with designed properties. This opens a similar 
degree of freedom to materials design as superlattices did for PCM’s, only with a wider 
range of “building blocks”. 

Designing these materials will be a challenge as the material property dependence on 
the layer sequence is likely to be very complex. Genetic algorithms (GA) have been 
used successfully in various tasks ranging from protein folding[T5] to crystal lattice 
optimization[TTj. GA’s generally work well in situations where the searched energy 
landscape has many minima that are separated by significant energy barriers and/or 
complicated reaction paths. This is because the stochastic nature of GA’s means they 
are less likely to converge to the local minima than deterministic methods [l6]. The ran¬ 
dom mutations and recombinations are the heart of GA’s. A search starts from a random 
sample of candidate solutions to the problem, their fitness to the solution is evaluated, 
and the properties of the best candidates are combined and mutated to create the next 
generation of candidates. 

Layered vdW heterostructures with Ge 2 Sb 2 Te 5 composition have show to perform 
higher phase-change efficiency than the alloy material of the same composition [12]. It is 
possible that they can be optimised even further. The aim of this paper is to show that 
GA-led design of 2D heterostructures can be used to hnd practical structures without 
the need to survey all the possible structures experimentally. The method is gener¬ 
ally applicable to many different problems, and as proof-of-concept it is applied to the 
optimization of phase-change superlattices. The optimization parameter used is the 
structure energy since it is relatively simple to calculate, and can be used to assess the 
stability of the structures grown. 

The tellurium sticking coefficient decreases at temperatures above 350 °G[2|. This 
makes it the highest practical deposition temperature for tellurium superlattices. There¬ 
fore a structure is estimated to be sufficiently stable for practical purposes if its energy 
is less than ksT = 53.7 meV/at. above the lowest structure energy found. This will 
yield a map of potential structures which can then be refined to find a set of desired 
properties such as optical contrast or resistivity change between states. As different 
vdW-heterostructures have a wide range of potential applications|l|, we envisage this 
method impacting on a range of fields from photonics to high-Tc superconductors |9|. 

2 Computational methods 

The genetic algorithm procedure is shown in figure the first generation of superlattice 
candidates was randomized, and subsequent generation candidates were combined or 
mutated from the previous generation. The algorithm uses permutations of 1-Natoms 
as intrinsic variables. This allows the adaptation of the crossover operator (OX) used 
in travelling salesman problems, which preserves the relative order of the variables|3| 
while combining genes from parent candidates to a child candidate. This speeds up 
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Figure 1: Flowchart of the algorithm 


the algorithm since a the order of the atomic layers determines the fitness so preserv¬ 
ing good partial sequences has a tendency to produce child candidates with a good 
fitness. The permutations are mapped to Ge, Sb and Te atomic sequences according 
to the desired material composition. The program starts by generating a population 
of permutations, for which it then calculates the fitness function with Vienna Ab initio 
Simulation Package (VASP)[6]. After calculations it reads the htness values, sorts the 
population accordingly, and creates the next generation of candidates. The VASP jobs 
can be executed sequentially, or submitted to a load leveler. In both cases it is possible 
to run the algorithm automatically up to a specified number of iterations. 

The present results were obtained with a population of 20 candidates, each subse¬ 
quent generation consisted of two elite candidates carried over from the previous gen¬ 
eration “as is” to preserve the best results throughout the run, 10 candidates produced 
from previous generation by crossover, and 8 candidates produced with mutation. The 
possible parent candidates for the next generation were the top 10 candidates. 

We used the final energy in the VASP geometry optimization with hxed X and 
Y coordinates of the atoms, and fixed simulation cell as the fitness function. Essen¬ 
tially the geometry optimization optimized the spacings between the atomic layers while 
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making sure that the order did not change. The simulation cell was hexagonal with 
a = b = 4.25 A, and c = Natoms x 1.9372 A, which resulted in the experimental crystal 
phase density of Ge 2 Sb 2 Te 5 . The initial atomic positions were laterally in the fractional 
coordinates (0,0), (2/3,1/3), and (1/3, 2/3), where the atoms were placed cyclically. The 
simulations were performed with 240 eV plane-wave cutoff energy, PBEsolfTO] exchange- 
correlation functional, PAW pseudopotentials [7], periodic boundary conditions, and an 
8x8x2 P centered Monkhorst-Pack k-point grid. 

An ordered crossover 1 (0X1) algorithm was used, which is good at preserving the 
sequence of genes. The parents, pi and p2, and two cutoff points were selected randomly. 
The genes of pi from first cutoff point to the second cutoff point were transfered to the 
offspring, and the remaining values were taken sequentially from the p2 starting from 
the gene immediately after cutoff point 2, and omitting the values already present. The 
mutation algorithm selected a random parent, and swapped two of its genes at random 
to produce a mutated offspring. 

When generating the initial population or the subsequent generations, it was required 
that the permutation for each candidate resulted in a unique atomic sequence for that 
generation. For example, if the mutation algorithm would pick an elite candidate, which 
is carried over to next generation unchanged, and swap two Te-atoms around in it, the 
resulting permutation would be rejected since it would result in an identical structure 
as a permutation already present in that generation. Cyclic permutations and reversed 
sequences were also rejected, since fitness function calculation used periodic boundary 
conditions, which deems such structures are identical. Whenever a generated candidate 
was rejected based on this condition, the generation procedure (randomization, com¬ 
bination or mutation) was repeated until a candidate unique for that generation was 
generated. 

3 Results 

Figure [^a) shows the energies of all candidates for each iteration. The energy of the 
lowest structure decreases by 0.8 eV over the course of the run, and reaches its lowest 
value after 42 iterations. The variety of structural energies for each iterations shows 
that the algorithm is trying structures away from the current minima in order to find 
better structures. The lowest energy structures are visualized in figure [^b), the common 
features of these are strong A-B-A-B alternation (A: Ge, Sb; B: Te), and separation of 
GeTe and Sb2Te3. The unavoidable Te-Te vdW interfaces favor Sb neighbours over Ge 
neighbours. 

A structure with the lowest energy is found, and it has an unusual atomic sequence 
that resembles two layers of quintuple-Sb 2 Te 3 stacked with one and three bilayers of 
GeTe. This means that the lowest energy structure cannot be modelled with 9 atom 
unit cell. We were able to grow this sequence onto a Si-substrate, and measure the 
XRD spectra. Figure]^ shows the comparison between measured XRD spectra, and the 
spectra calculated from the simulation coordinated using Materials Studio (version 8.0) 
of Accelrys Inc. The agreement is good, most peaks are present in both spectra, and 


4 


4 


3.5 
3 -! 

25 il I • • ■ . • 

.1* !;»•* . I . •. ••.«• ••! !* i 
:il . I'!" 

15 ! 

I •• •••> * It**! * 

*!tl*i*i**^itlii* ****t*******!: 

' L::;; 

C. •« • . 

•• ••• 






. .::t., 


5 10 15 20 25 30 35 40 45 50 

Generation # 

(a) Candidate Energies 



(b) Lowest energy candidate structures 


Figure 2: (a) Energies of all candidates as a function of the generation number shown in 
blue, lowest energy evolution is hilighted in red. (b) Visualizations of the lowest energy 
structures for each iteration. Colors: germanium in tan, antimony in yellow, tellurium 
in cyan. 


only one small peak at ~29° is present only in the measurement, and one peak at ~39° 
is present only in the simulated spectra. In addition to the lowest energy structure we 
obtained a large number of relatively low energy structures that are therefore practi¬ 
cal, and can be investigated further for phase-change, atomic switching, photonics, and 
memory applications. To elaborate on this, we looked for the lowest energy structure 
that contains Ge-Te-Te-Ge sequence, as that is regarded as the location where the Ge 
atomic switching is confined. This structure is referred to as “switchable” structure, and 
it is visualized in figure]^ This structure is 8.52 meV/at. higher than the lowest energy 
structure. 

We compared our lowest energy structure with some of the Ge 2 Sb 2 Te 5 crystal and 
superlattice structures reported in literature [T3l [5l [H] by relaxing the respective struc¬ 
tures and their simulation cells. Our results are given in table and show that the 
structures found by GA are very close to the energetically most stable stucture, which is 
that suggested by Kooi and Hosson[5]. The switchable structure is higher in energy by 
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Figure 3: Calculated XRD spectra for the lowest energy superlattice structure (blue), 
and the measured XRD spectra for the sputtered sample (red) 


Table 1: Energies and simulation cell parameters of different layer orders of Ge 2 Sb 2 Te 5 . 
The overall lowest energy sequence, and the lowest 


Kooi[5] 

E (meV/at.) 
0.00 

a (A) 
4.30 

c(A) 

34.85 

Present work 

0.92 

4.30 

35.07 

Ferro GeTe[T3] 

8.12 

4.28 

35.34 

Present work (switchable) 

9.44 

4.28 

35.36 

Inverse Petrov [13] 

13.14 

4.21 

38.23 

Petrov [TT] 

17.41 

4.27 

35.66 


9.44 meV/at, and is energetically close to the ferro-GeTe structure. It should be noted 
that the all but inverse Petrov structures have similar lattice constants, and could likely 
coexist in an imperfect crystal. 

4 Conclusions 

In this paper, a genetic algorithm to efficiently design layer sequences for vdW het¬ 
erostructures was introduced and applied to Ge 2 Sb 2 Te 5 superlattices. The algorithm 
found structures with energies comparable to the energetically most favorable known 
Ge 2 Sb 2 Te 5 ground state structures along with “switchable” structures which contain 
Ge-Te-Te-Ge layer sequence. The best structure was grown on a silicon substrate and 
its XRD pattern agreed well with the spectra computed for the modeled structure. This 
lends further support that the structures generated by the algorithm are practical. 

The fitness function in the algorithm can in practice be any property which can be 
computed for a structure. This makes it versatile in finding structures or structure com- 
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Figure 4: The structure of the lowest energy structure with Ge-Te-Te-Ge sequence found. 
Golors as in fig. j^b) 


binations that can be used in various applications ranging from photonics to memory 
devices. The algorithm can be used as a too to help to identify structure-property corre¬ 
lation or structural motifs that can be used to design new van-der-Waals heterostructure 
superlattices with interesting properties. 
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